SOME ASPECTS OF LARGE TIME BEHAVIOR OF 
THE HEAT KERNEL: AN OVERVIEW WITH 
PERSPECTIVES 



YEHUDA PINCHOVER 

_ . Dedicated to Professor Michael Demuth on the occasion of his 65th birthday 

<D 

| Abstract. We discuss a variety of developments in the study 

of large time behavior of the positive minimal heat kernel of a 
time independent (not necessarily symmetric) second-order para- 
bolic operator defined on a domain M C Mr, or more generally, on 
a noncompact Ricmannian manifold M . Our attention is mainly 
focused on general results in general settings. 
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1. Introduction 

The large time behavior of the heat kernel of a second-order parabolic 
operator has been extensively studied over the recent decades (see for 
example the following monographs and survey articles [H [TTJ [TTJ [18j 
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EQl [251 [291 [371 EH EHJ [571 [591 E21 [631 El] , and references therein) . The 
purpose of the present paper is to review a variety of developments 
in this area, and to point out a number of their consequences. Our 
attention is mainly focused on general results in general settings. Still, 
the selection of topics in this survey is incomplete, and is according 
to the author's working experience and taste. The reference list is far 
from being complete and serves only this expose. 

Let P be a general linear, second-order, elliptic operator defined on 
a domain M C M d or, more generally, on a noncompact, connected, 
Riemannian manifold M of dimension d > 1. Denote the cone of 
all positive solutions of the equation Pu = in M by Cp{M). The 
generalized principal eigenvalue is defined by 

A = A (P,M) := sup{A G R | C P _ A (M) ^ 0}. 

Throughout this paper we always assume that Ao > — oo. 

Suppose that Ao > 0, and consider the (time- independent) parabolic 
operator 

(1.1) Lu := d t u + P(x,d x )u (x, t) G M x (0, oo). 

We denote by %p(M x (a, b)) the cone of all nonnegative solutions of 
the parabolic equation 

(1.2) Lu = mMx(a,b). 

Let kp(x,y,t) be the positive minimal heat kernel of the parabolic 
operator L on the manifold M. By definition, for a fixed y G M, the 
function (x,t) >->• kp(x,y,t) is the minimal positive solution of the 
equation 

(1.3) Lu = inMx(0,oo), 

subject to the initial data S y , the Dirac distribution at y G M. It can 
be easily checked that for A < A , the heat kernel kp_ x of the operator 
P — A on M satisfies the identity 

(1.4) k^_ x (x,y,t) = e xt k^(x,y,t). 

So, it is enough to study the large time behavior of kp_ Xo , and therefore, 
in most cases we assume that \q(P, M) = 0. Note that the heat kernel 
of the operator P*, the formal adjoint of the operator P on M, satisfies 
the relation 

k£(x,y,t) = k¥(y,x,t). 

Here we should mention that many authors derived upper and lower 
Gaussian bounds for heat kernels of elliptic operators on M := M. d , 
or more generally, on noncompact Riemannian manifolds M. As a 
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prototype result, let us recall the following classical result of Aronson 
®- 

Example 1.1. Let P be a second-order uniformly elliptic operator 
in divergence form on IR d with real coefficients satisfying some gen- 
eral boundedness assumptions. Then the following Gaussian estimates 
hold: 

-d/2 p ^ (_ r J x -y\ z _,, U A < 



C 1 (4 7 rt)- c( / i exp l-C 2 ^ t yi -u x t\ < kf(x,y,t) 

< C7 3 (47rt)- d/2 exp ^-cJ-^—^- + co 2 t S j V{x,y,t) ER d xR d xR + . 

However, since Gaussian estimates of the above type in general are 
not tight as t — > oo, such bounds do not provide us with the exact 
large time behavior of the heat kernel, let alone strong ratio limits of 
two heat kernels. 

In spite of this, and as a first and rough result concerning the large 
time behavior of the heat kernel, we have the following explicit and 
useful formula 

(1.5) lim lQ g fc P = _A . 

<->oo t 

We note that (11.51) holds in the general case, and characterizes the 
generalized principal eigenvalue Ao in terms of the large time behavior 
of log kp(x, y, t). The above formula is well known in the symmetric 
case, see for example [31] and [27] Theorem 10.24]. For the proof in 
the general case, see Corollary 16.11 

To get a more precise result one should introduce the notion of criti- 
cality. We say that the operator P is subcritical (respectively, critical) 
in M if for some x ^ y, and therefore for any x ^ y, x, y G M, we have 
(1.6) 

/oo / roc 

kp(x, y, r) dr < oo I respectively, / kp(x, y, r) dr 

It follows from the above definition that, roughly speaking, the heat 
kernel of a subcritical operator in M "decays" faster as t — > oo than 
the heat kernel of a critical operator in M. This rule of thumb will be 
discussed in Section [HI 

If P is subcritical in M, then the function 



oo 



poo 

;i.7) G^(x,y):= / k?(x,y,T)dr 

Jo 
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is called the positive minimal Green function of the operator P in M. 

It follows from (jl.5p that for A < Ao, the operator P — A is subcritical 
in M. Clearly, P is subcritical (respectively, critical) in M, if and only 
if P*, is subcritical (respectively, critical) in M. Furthermore, it is well 
known that if P is critical in M, then Cp(M) is a one-dimensional cone, 
and any positive supersolution of the equation Pu = in M is in fact a 
solution. In this case, the unique positive solution ip G Cp(M) is called 
Agmon ground state (or in short ground state) of the operator P in M 
[U SHI \M\- We denote the ground state of P* by ip*. 

The following example demonstrates two prototype behaviors as t — > 
oo of heat kernels corresponding to two particular classical cases. 

Example 1.2. I. Let P = -A, M = R d , where d > 1. It is well 

known that Ao(— A,R d ) = 0, and that the heat kernel is given by the 
Gaussian kernel. Hence, 

(x,y,t) = k^ d A (x,y,t) 

12 



exp 



\x y | \ dii c\ 

- M i — : I ~ t ' — > 0. 

(47rt)d/2 ^ ^ 4t y t^oo t^oo 

So, the rate of the decay depends on the dimension, and clearly, —A is 
critical in IR d if and only if d — 1, 2. Nevertheless, the above limit is 
in any dimension. 

2. Suppose that P is a symmetric nonnegative elliptic operator with 
real and smooth coefficients which is defined on a smooth bounded do- 
main M (s WL d , and let {<Pn}^Lo be the complete orthonormal sequence 
of the (Dirichlet) eigenf unctions of P with the corresponding nonde- 
creasing sequence of eigenvalues {Anj-^Q. Then the heat kernel has the 
eigenfunction expansion 

oo 

(1.8) kp*(x,y,t) = J2 e ~ Xnt( Pn(x) Vn (y). 

n=0 

Hence, 

oo 

e Xot k^{x,y,t) = Ve^-^Vn^KO/) -> Mx)My) > 0. 

< * t— s-no 



n=0 



Therefore, the operator P— Ao is critical in M, t 1 log kp_ Xo (x, y, t) 
0, but kp_ Xo does not decay as t — > oo. 



t— >oo 
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3. Several explicit formulas of heat kernels of certain classical oper- 
ators are included in [TUl EHl EZ] , each of which either tends to zero or 
converges to a positive function as t — > oo. 

The characterization of Ao in terms of the large time behavior of the 
heat kernel, given by (II. 5p . provides us with the asymptotic behavior 
of \ogkp as t — > oo but not of kp itself. Moreover, ( 11. 5ft does not 
distinguish between critical and subcritical operators. In the first part 
of the present article we provide a complete proof that 

lim e Xot k^{x,y,t) 

t—>oo 

always exists. This basic result has been proved in two parts in jH] 
and |51j . and here, for the first time, we give a comprehensive and a bit 
simplified proof (see also [5l[12l[28l[^[32l[Ml[ml[62l[66] and references 
therein for previous and related results). We have: 

Theorem 1.3 ( |^8l IBTj ). Let P be an elliptic operator defined on M, 
and assume that A (-P, M) > 0. 

(i) : The subcritical case: If P — X is subcritical in M, then 

lim e Xot k^(x,y,t) = 0. 

t— >oo 

(ii) : The positive-critical case: If P — X critical is in M, and 
the ground states and tp* of P — Ao and P* — Ao, respectively, 
satisfy ip*ip G L l (M), then 

lime^( X ,y,t) = ^ X)V ^ y) 

t—>oo 



(Hi): The null-critical case: If P — Ao is critical in M, and 
the ground states <p and <p* of P — Ao and P* — X , respectively, 
satisfy tp*tp £ L l (M), then 

lim e Xot k^{x,y,t) = 0. 



t— >oo 



Moreover, for X < Ao, let G P _ x (x,y) be the minimal positive Green 
function of the elliptic operator P — X on M. Then the following 
Abelian-Tauberian relation holds 

(1.9) lim e Xot k^{x, y, t) = lim (A - X)G™(x, y). 

t— s-oo A/'Ao 

The outline of the present paper is as follows. In Section |2] we pro- 
vide a short review of the theory of positive solutions and the basic 
properties of the heat kernel. The proof of Theorem 11.31 is postponed 
to Section [5] since it needs some preparation. It turns out that the 
proof of the null-critical case (given in [51 J) is the most subtle part of 
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the proof of Theorem 11.31 It relies on the large time behaviors of the 
parabolic capacitory potential and of the heat content that are studied 
in Section [3j and on an extension of Varadhan's lemma that is proved 
in Section H] (see Lemma H~T!) . Section [6] is devoted to some applications 
to Theorem 11.31 In particular, Corollary 16.71 seems to be new. 

The next two sections deal with ratio limits. In Section [7J we discuss 
the existence of the strong ratio limit 

k?(x,y,t) 



t^oo k^(x ,x ,t) 

(Davies' conjecture), and in Section |8] we deal with the conjecture that 
(1.10) lim^ -( = Vx,yeM, 

where P + and P are respectively subcritical and critical operators in 
M. Finally, in Section [9] we discuss the equivalence of heat kernels of 
two subcritical elliptic operators (with the same Ao) that agree outside a 
compact set. Unlike the analogous question concerning the equivalence 
of Green functions, the problem concerning the equivalence of heat 
kernels is still "terra incognita" . Solving this problem will in particular 
enable us to study the stability of the large time behavior of the heat 
kernel under perturbations. 

The survey is an expanded version of a talk given by the author at 
the conference "Mathematical Physics, Spectral Theory and Stochastic 
Analysis" held in Goslar, Germany, September 11-16, 2011, in honor 
of Professor Michael Demuth. We note that most of the results in 
sections [3J-E1 originally appeared in [JEHS], those of Section [7] appeared 
in [52], while those in sections [HHS1 originally appeared in 



2. Preliminaries 

In this section we recall basic definitions and facts concerning the 
theory of nonnegative solutions of second-order linear elliptic and par- 
abolic operators (for more details and proofs, see for example [ 5 4j ) . 

Let P be a linear, second-order, elliptic operator defined in a non- 
compact, connected, C 3 -smooth Riemannian manifold M of dimension 
d. Here P is an elliptic operator with real and Holder continuous coef- 
ficients which in any coordinate system (U; x±, . . . , x&) has the form 

d d 

(2.1) P(x, d x ) = - ^ a-ij(x)didj + ^ H x ) d i + c(x), 

i,j=l i=l 
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where di = d/dxi. We assume that for every x G M the real quadratic 
form 

d 

(2.2) ">M)^.r Z = (Zi,...,Zd)eR d 

is symmetric and positive definite. 

Throughout the paper we always assume that Aq(P, M) > 0. In this 
case we simply say that the operator P is nonnegative in M, and write 
P > in M. So, P > in M if and only if the equation Pu = admits 
a global positive solution in M. We note that in the symmetric case, 
by the Agmon-Allegretto-Piepenbrink theory [I], the nonnegativity of 
P in M is equivalent to the nonnegativity of the associated quadratic 
form on C °°(M). 

We consider the parabolic operator L 

(2.3) Lu = u t + Pu onMx(0,oo), 

and the corresponding homogeneous equation Lu = in M x (0, oo). 

Remark 2.1. We confine ourselves to classical solutions. However, 
since the results and the proofs throughout the paper rely on standard 
elliptic and parabolic regularity, and basic properties and results of po- 
tential theory, the results of the paper are also valid for weak solutions 
in the case where the elliptic operator P is in divergence form 

d 

<iij(x)djU + ubi(x) 
i=i Li=i J j=i 

with coefficients which satisfy standard local regularity assumptions 
(as for example in Section 1.1 of [B]). The results are also valid in the 
framework of strong solutions, where the strictly elliptic operator P is 
of the form (12.11) and has locally bounded coefficients; the proofs differ 
only in minor details from the proofs given here. 

We write Q\ d ^2 if ^2 is open, is compact and £1% C ^2- Let 
f,g€ C(Q) be nonnegative functions; we use the notation / x g on 
if there exists a positive constant C such that 

C^g(x) < f(x) < Cg(x) for all x e ft. 

Let {Mj}'jL 1 be an exhaustion of M, i.e. a sequence of smooth, 
relatively compact domains in M such that Mi 7^ 0, Mj <s Mj + i and 
\Jf =l Mj = M. For every j > 1, we denote M* = M\c\{Mj). Let = 
M U {00} be the one-point compactification of M. By a neighborhood 
of infinity in M we mean a neighborhood of 00 in M^, that is, a set of 
the form M\K, where K is compact in M. 



(2.4) P u = -J2di 



+ bi{x)diU + c(x)u, 
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Assume that P > in M. For every j > 1, consider the Dirichlet 
heat kernel kp ] (x,y,t) of the parabolic operator L — d t + P in Mj. 
So, for every continuous function / with a compact support in M, the 
function 

u(x,t) := kp j (x,y,t)f(y)dy 
Jm 3 

solves the initial-Dirichlet boundary value problems 

Lu = in Mj x (0, oo), 

(2.5) u = on dMj x (0,oo), 

u = f on Mj x {0}. 

By the generalized maximum principle, {kp J (x, y, t)}'jL 1 is an increas- 
ing sequence which converges to kp(x,y,t), the positive minimal heat 
kernel of the parabolic operator L in M. 

The main properties of the positive minimal heat kernel are summa- 
rized in the following lemma. 

Lemma 2.2. Assume that P > in M. The heat kernel kp(x,y,t) 
satisfies the following properties. 

(1) Positivity: kp(x, y, t) > for allt>0 and x,y G M. 

(2) For any \ < \ we have 

(2.6) kp*_ x (x,y,t) = e xt kp*(x,y,t). 

(3) The heat kernel satisfies the semigroup identity 
(2.7) 

/ kp(x, z, t)kp(z, y, r) dz = kp(x, y,t + r) Vt, r > 0, and x,y e M. 
Jn 

(4) Monotonicity: If M 1 C M, and Vi > V^, £/ien 

(5) For any fixed y G M (respectively x <E M), kp(x,y,t) solves the 
equation u t + P(x, <9 x )-u = (respectively u t + P*(y, c^w = 

in M x (0, oo). 

(6) S&ew; product operators: Let M = Mi x M2 fre a product of two 
manifolds, and denote x = (x±,X2) G M = Mi x M 2 . Consider 
a skew product elliptic operator of the form P := Pi<8>l2+Ii<8>P2, 
where Pi is a second-order elliptic operator on Mi satisfying the 
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assumptions of the present paper, and Li is the identity map on 
Mi, i = 1,2. Then 

(2.8) k¥(x,y,t) = k^(x l ,y l ,t)k^(x 2 ,y 2 ,t). 

(7) Eigenfunction expansion: Suppose that P is a symmetric ellip- 
tic operator with (up to the boundary) smooth coefficients which 
is defined on a smooth bounded domain M. Let {(p n }^ =0 be the 
complete orthonormal sequence of the (Dirichlet) eigenf unctions 
of P with the corresponding nondecreasing sequence of eigenval- 
ues {Anj^Q. Then the heat kernel of P in M has the eigen- 
function expansion 

oo 

(2.9) k™(x, y,t) = J2 e- Xnt M*)My)- 

n=0 

(8) Let v E C P {M) and v* E C P *{M). Then we have [Ml S9] 

Jm Jm 

Moreover, (by the maximum principle) either 

/ k¥(x,y,t)v(y)dy < v(x) V(x, t) E M x (0, oo), 
Jm 

or 

/ kp(x,y,t)v(y)dy = v(x) V(x, t) E M x (0, oo), 
Jm 

and in the latter case v is called an invariant solution of the 
operator P on M (see for example [211 ESI SHI IB]/ 

Remark 2.3. 1. If there exists v E Cp(M) such that v is not invariant, 
then the positive Cauchy problem 

Lu = 0,u>0 onMx(0,oo), u(x,0) = xEM 

does not admit a unique solution. 

2. An invariant solution is sometimes called complete. If the constant 
function is an invariant solution with respect to the heat operator, 
then one says that the heat operator conserves probability, and the 
corresponding diffusion process is said to be stochastically complete 
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Definition 2.4. Let to be a positive solution of the equation Pu = 
in M \ K, where K d= M. We say that w is a positive solution of the 
equation Pu = of minimal growth in a neighborhood of infinity in 
M if for any v G C(cl(M*)) which is a positive supersolution of the 
equation Pu = in M* for some j > 1 large enough, and satisfies 
w <v on dM*, we have w < v in M* [Tlli8|[54]. 

In the subcritical (respectively, critical) case, the Green function 
Gp(-,y) (respectively, the ground state <p) is a positive solution of the 
equation Pu = of minimal growth in a neighborhood of infinity in 
M. Recall that if A < Ao, then P — A is subcritical in M. In particular, 
if P is critical in M, then Ao = 0. 

Next, we recall the parabolic Harnack inequality. We denote by 
Q(x , t , R, t) the parabolic box 

Q(x , t , R, 9) := {(x, t) 6 M x R | p(x, x ) < R,t e (to, t + ^ 2 )}, 

where p is the given Riemannian metric on M. We have: 

Lemma 2.5 (Harnack inequality). Let u be a nonnegative solution of 
the equation Lu = in Q(xq, to, R, 9), and assume that 9 > 1 and 
0<R<Ro. Then 

(2.11) w(x ,to + J R 2 ) < Cu(x,t + 9R 2 ) 

for every p(x, Xo) < R/2, where C — C(L, d, Ro, 9). Moreover, C varies 
within bounded bounds for all 9 > 1 such that < e < (9 — 1) _1 < M. 

Let f G Cp(M) and G Cp*(M). Then the parabolic Harnack 
inequality (12. lip , and (I2.10p imply the following important estimate 

(2.12) kg{x,y,t) < c^Mx), and k^{x,y,t) < c 2 (x)v*(y) 

for all x, y G M and t > 1 (see [48, (3.29-30)]). Recall that in the crit- 
ical case, by uniqueness, v and v* are (up to a multiplicative constant) 
the ground states p and (p* of P and P* respectively. Moreover, it is 
known that the ground state <p (respectively, <p*) is a positive invariant 
solution of the operator P (respectively, P*) in M. So, 

We distinguish between two types of criticality. 

Definition 2.6. A critical operator P is said to be positive- critical in 
M if <p> G L X (M), and null-critical in M if v?> ^ L X (M). 
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Remark 2.7. Let 1 be the constant function on M, taking at any 
point x G M the value 1. Suppose that PI = 0. Then P is sub- 
critical (respectively, positive-critical, null-critical) in M if and only if 
the corresponding diffusion process is transient (respectively, positive- 
recurrent, null-recurrent). For a thorough discussion of the probabilis- 
tic interpretation of criticality theory, see [54J. 

In fact, in the critical case it is natural to use the well known (Doob) 
/i-transform with h = <p, where ip is the ground state of P. So, 

P^u := —P(ifu) and therefore kp lfi (x,y,t) = r kp (x,y,t)<p(y). 



Clearly, P v is an elliptic operator which satisfies all our assumptions. 
Note that P v is null-critical (respectively, positive-critical) if and only if 
P is null-critical (respectively, positive-critical), and the ground states 
of P v and (P 1 ^)* are 1 and <p*ip, respectively. Moreover, 

lim k PV (x, y, t) — if and only if lim k P (x, y, t) = 0. 

t— >oo t— >oo 

Therefore, in the critical case, we may assume that 

(A) PI = 0, and P is a critical operator in M. 

It is well known that on a general noncompact manifold M, the 
solution of the Cauchy problem for the parabolic equation Lu = is not 
uniquely determined (see for example [30] and the references therein). 
On the other hand, under Assumption (A), there is a unique minimal 
solution of the Cauchy problem and of certain initial-boundary value 
problems for bounded initial and boundary conditions. More precisely, 

Definition 2.8. Assume that PI = 0. Let / be a bounded continuous 
function on M. By the minimal solution u of the Cauchy problem 



(p{x) 



Lu 



u 



in M x (0, oo) 
/ on M x {0}, 



we mean the function 



(2.14) 




Note that fl2TT2l implies that u in fF2TT4l is well defined. 



Definition 2.9. Assume that PI = 0. Let B d Mi be a smooth 
bounded domain such that B* := M \ c\(B) is connected. Assume 
that / is a bounded continuous function on B*, and g is a bounded 
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continuous function on OB x (0, oo). By the minimal solution u of the 
initial-boundary value problem 

Lu = in B*x (0,oo), 

(2.15) u = g on OB x (0,oo), 

u = f on B* x {0}, 

we mean the limit of the solutions Uj of the following initial-boundary 
value problems 

Lu = in (S* nMj) x (0,oo), 
u = g on dB x (0, oo), 
u = on <9Mj x (0,oo), 
u = f on (5* n M,-) x {0}. 

Remark 2.10. It can be easily checked that the sequence {uj} is in- 
deed a converging sequence which converges to a solution of the initial- 
boundary value problem (I2.15p . 

Next, we recall some results concerning the theory of positive solu- 
tions of elliptic equations that we shall need in the sequel. 

The first result is a Liouville comparison theorem in the symmetric 
case. 

Theorem 2.11 ([53]). Let Po and P\ be two symmetric operators de- 
fined on M of the form 

(2.16) PjU = -m^divimjAjVu) + Vju j = 0, 1. 

Assume that the following assumptions hold true. 

(i) The operator Po is critical in M. Denote by (p G Cp (M) its 
ground state. 

(ii) Pi > in M, and there exists a real function ip 6 Hl oc (M) such 
that tp + 7^ 0, and P^ < in M , where u + (x) := max{0, u(x)}. 

(iii) The following matrix inequality holds 
(2.17) 

(ip+) 2 (x)mi(x)Ai(x) < Cip 2 (x)mo(x)A (x) for a.e. x E M, 

where C > is a positive constant. 

Then the operator Pi is critical in M, and ip is its ground state. In 
particular, dimCp x (M) = 1 and X (Pi,M) = 0. 

In the sequel we shall also need to use results concerning small and 
semismall perturbations of a subcritical elliptic operator. These notions 
were introduced in [26] and [23] respectively, and are closely related to 
the stability of Cp(fl) under perturbation by a potential V. 
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Definition 2.12. Let P be a subcritical operator in M, and let V be 
a real valued potential defined on M. 

(i) We say that V is a small perturbation of P in M if 

(2.18) ito/sup / gycy^^Lo, 

(m,) V is a semismall perturbation of P in M if for some Xo € M we 
have 

(2.19) Km { sup / = o. 
' ' s {:, M] Jai; Gp{x ,y) J 

Recall that small perturbations are semismall [33]. For semismall 
perturbations we have 

Theorem 2.13 ([321 1351 HI])- Let P be a subcritical operator in M. 
Assume that V = V + — V- is a semismall perturbation of P* in M 
satisfying V- ^ ; where V±(x) = max{0, ±V(x)}. 

Then there exists ao > such that P a := P + aV is subcritical in M 
for all < a < «o and critical for a = ao. 

Moreover, let <p be the ground state of P + aoV and let yo be a fixed 
reference point in Mi . Then for any < a < ao 

F^( Jo ) in Ml 

where the equivalence constant depends on a. 

3. CAPACITORY POTENTIAL AND HEAT CONTENT 

Our first result concerning the large time behavior of positive solu- 
tions is given by the following simple lemma that does not distinguish 
between null-critical and positive-critical operators. 

Lemma 3.1. Assume that PI = and that P is critical in M. Let 
B := B(xq, 5) CC M be the ball of radius 5 centered at xq, and suppose 
that B* = M \ cl(B) is connected. Let w be the heat content of B* , 
i.e. the minimal nonnegative solution of the following initial-boundary 
value problem 

Lu = in B* x (0, oo), 

(3.1) u = ondB x (0,oo), 

u = 1 on B* x {0}. 

Then w is a decreasing function oft, and lim^oo w(x,t) = locally 
uniformly in B* . 
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Proof. Clearly, 

(3.2) w(x,t)= f kf(x,y,t)dy<[ kp*{x, y, t) dy = 1. 

J B* J M 

It follows that < w < 1 in B* x (0, oo). Let e > 0. By the semigroup 
identity and ( 13. 2p . 

w(x,t + e) = / kp (x,y,t + e) dy = 
Jb* 

(3.3) ^ k B p* (x, z, t)k B p* {z, y, e) dz^j dy = 

(3.4) kf (x, z, t) Qf^ kf (z, y, e) dyj dz < 

/ kp (x, z, t) dz = w(x, t). 
Jb* 

Hence, if is a decreasing function of t, and therefore, lim^oo w(x, t) ex- 
ists. We denote v(x) := ]im. t _ >00 w(x, t). Note that the above argument 
shows that even in the subcritical case w is a decreasing function of t. 

For t > 0, consider the function u(x, i; r) := w(x,t + r), where t > 
—t. Then ^(x, t; r) is a nonnegative solution of the parabolic equation 
Lu = in B* x (— r, oo) which satisfies u = on 95 x (— r, oo). By a 
standard parabolic argument as r — > oo, any converging subsequence 
of this set of solutions converges locally uniformly to a solution of the 
parabolic equation Lu = in B* x K which satisfies m = on 95 x IR. 
Since linv-Kx, u(x, t; r) = lim r ^.oo w(x, r) = v(x), the limit does not 
depend on t, and v is a solution of the elliptic equation Pu = in B*, 
and satisfies v = on OB. Furthermore, < v < 1. 

Therefore, 1 — t> is a positive solution of the equation Pu = in 
5* which satisfies n = 1 on <9-B. On the other hand, it follows from 
the criticality assumption that 1 is the minimal positive solution of the 
equation Pu = in B* which satisfies u = 1 on OB. Thus, 1 < 1 — v, 
and therefore, v = 0. □ 



Definition 3.2. Let B := B(x ,5) CC M. Suppose that B* 
M \ c\{B) is connected. The nonnegative (minimal) solution 



v(x,t) = 1-1 kp*(x,y,t)dy 

B* 
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is called the parabolic capacitory potential of B*. Note that v is indeed 
the minimal nonnegative solution of the initial-boundary value problem 

Lu = in B* x (0,oo), 

(3.5) u = 1 on^Bx (0,oo), 

u = on 5* x {0}. 

Corollary 3.3. Under the assumptions of Lemma \3.1\ the parabolic 
capacitory potential v of B* is an increasing function oft, and we have 
lim^oo v(x, t) = 1 locally uniformly in B* . 

Proof. Using an exhaustion argument it is easily verified that 

(3.6) v(x, t) = 1 — / kp*{x, y, t) dy = 1 — w(x, t), 

J B* 

where w is the heat content of B*. Therefore, the corollary follows 
directly from Lemma [3.11 □ 

4. VARADHAN'S LEMMA 

Varadhan's celebrated lemma (see, [621 Lemma 9, p. 259] or [54, 
pp. 192-193]) deals with the large time behavior of minimal solutions of 
the Cauchy problem with bounded initial data (assuming that PI = 0). 
It turns out that the limit as t — > oo of such solutions might not exist, 
but, under further conditions, the spacial oscillation of the solution 
tends to zero as t — > oo. In the present section we slightly extend 
Varadhan's lemma (Lemma 14. ip . This extended version of the lemma 
is crucially used in the proof of the null-critical case in Theorem 11.31 

Varadhan proved his lemma for positive-critical operators on M. d us- 
ing a purely probabilistic approach (ours is purely analytic). Our key 
observation is that the assertion of Varadhan's lemma is valid in our 
general setting under the weaker assumption that the skew product 
operator P :— P <g) / + / <S> P is critical in M := M x M, where / is 
the identity operator on M. Note that if P is critical in M, then P is 
critical in M. On the other hand, if P is positive-critical in M, then 
P is positive-critical in M. Moreover, if P is subcritical in M, then by 
part (i) of Theorem 11.31 the heat kernel of P on M tends to zero as 
t — > oo. Since the heat kernel of P is equal to the product of the heat 
kernels of its factors (see Lemma \2. 2 1) , it follows that if P is subcritical 
in M, then lim^oo kp(x, y, t) = 0. 

Consider the Riemannian product manifold M = M x M. A point in 
M is denoted by x = (x±, xq). By P Xi , i = 1, 2, we denote the operator 
P in the variable X{. So, P = P X1 + P X2 is in fact the above skew 
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product operator defined on M. We denote by L the corresponding 
parabolic operator. 

Lemma 4.1 (Varadhan's lemma [SI]). Assume that PI = 0. Suppose 
further that P is critical on M . Let f be a continuous bounded function 
on M , and let 

u(x,t)= f kp*(x,y,t)f(y)dy 

J M 

be the minimal solution of the Cauchy problem with initial data f on 
M. Fix K CC M. Then 

lim sup \u(xi, t) — u(x 2 , t)\ = 0. 
x!,x 2 eK 

Proof. Denote by u(x,t) := u(x\,t) — u(x 2 ,t). Recall that the heat 
kernel k(x,y,t) of the operator L on M satisfies 

(4.1) kf(x,y,t) = k^(x 1 ,y l ,t)k^ I (x 2} y 2 ,t). 
By (J2II3D and flUD, we have 

u(x, t) = u(xi,t) — u(x 2 , t) = 
/ kp (x u yi, t)f (yx) dyx - / kp {x 2 ,y 2 ,t)f(y 2 ) dy 2 = 

J M J M 

[ [ k K p I {x l ,y u t)k^ I {x 2) y 2 Mf{yi)-f{y2))dy l dy 2 = 

J M J M 

[ kf(x,y,t)(f( yi )-f(y 2 ))dy. 

J M 

Hence, u is the minimal solution of the Cauchy problem for the equa- 
tion Lw = on M x (0, oo) with the bounded initial data f(x) : = 
f{xi) — f(x 2 ), where x G M. 

Fix a compact set K CC M and Xq G M \ K, and let £ > 0. Let 
B := B((x ,x Q ),S) CC M \ K, where K = K x K, and 5 will be 
determined below. We may assume that B* = M\ cl(B) is connected. 
Then u is a minimal solution of the following initial-boundary value 
problem 

Lu = in B* x (0,oo), 

(4.2) u(x,t) = u(xi,t) — u(x 2 ,t) on dB x (0, oo), 
u(x,0) = f(xi)-f(x 2 ) on B* x {0}. 

We need to prove that lim^oo u{x,t) = 0. 

By the superposition principle (which obviously holds for minimal 
solutions), we have 

u(x, t) = ui(x, t) + u 2 (x, t) on B* x [1, oo), 
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where u± solves the initial-boundary value problem 

Lu\ = in B* x (1, oo), 

(4.3) Ui{x,t) = u(xi, t) — u(x 2 , t) on dB x (1, oo), 
Ul (x,0) = on B* x {1}, 

and u 2 solves the initial-boundary value problem 

Lu 2 = in B* x (1, oo), 

(4.4) u 2 (x,t) = on dB x (l,oo), 
u 2 (x,0) = u(xi, 1) — u(x 2 , 1) on B* x {1}. 

Clearly, \u{x,t)\ < 2|| /]|oo on M x (0,oo). Note that if x = (xi,x 2 ) G 
<9.B, then on M, distA/(a;i, £2) < 25. Using Schauder's parabolic interior 
estimates on M, it follows that if 5 is small enough, then 

\u(x, t)\ = \u(xi, t) — u(x 2 , t)\ < £ on dB x (1, 00). 

By comparison of u\ with the parabolic capacitory potential of B*, we 

obtain that 

(4.5) 

\v>i{%, t)\ < e — J kp*(x,y, t — 1) dy^j < e in £?* x (1, 00). 
On the other hand, 

(4.6) MM) I <2||/|| Q0 / kf(x,y,t-l)dy in B* x (1, 00). 

It follows from (14. 6p and Lemma 13.11 that there exists T > such that 

(4.7) \u 2 (x,t)\<e for all x G K and t > T. 

Combining (14. 5 p and (14.71) . we obtain that \u(xi,t) — u(x 2 ,t)\ < 2s for 
all x\, x 2 G K and t > T. Since e is arbitrary, the lemma is proved. □ 

5. Existence of \im t ^ 00 e x ° t k^(x,y,t) 

The present section is devoted to the proof of Theorem 11.31 which 
claims that lim^oo e Xot kp(x, y, t) always exists. Without loss of gen- 
erality, we may assume that Ao = 0. For convenience, we denote by k 
the heat kernel kp of the operator P in M. 

Proof of Theorem \1.3[ (i) The subcritical case: Suppose that P is sub- 
critical in Q. It means that for any x, y, G M we have 

/oo 
fc(x, y, t) dt < 00. 

Let x,y G Q be fixed, and suppose that k(x,y,t) does not converge 
to zero as t — > 00. Then there exist an increasing sequence tj — > 00, 
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tj+i — tj > 1, and £ > 0, such that k(x,y,tj) > e for all j > 1. Using 
the parabolic Harnack inequality ( 12. lip , we deduce that there exists 
C > (C may depend on x, y, Q, and P) such that 

(5.2) ! 3+ k{x,y,t)dt> Ce Vj > 2. 



But this contradicts (15.11) . 

(ii) The positive- critical case: Let P be a critical operator in M and 
let {tj}jt 1 C 1R be a sequence such that tj — > oo, and define: 

Uj(x,y,t) = k(x,y,t + tj). 

Fix x , ?/o e M. By (12~T2D we have 

(5.3) Uj(x,y ,t) < ci(y )(p(x), (uj(x ,y,t) < c 2 (x )^*(y). i 



for all x G M (y G M) and t E 1 and j > Using the parabolic 

Harnack inequality and a standard Parabolic regularity argument we 
may subtract a subsequence of {uj} (which we rename by {uj}) that 
converges locally uniformly to a nonnegative solution u(x,y,t) of the 
parabolic equations 

Lu = d t u + P(x, d x )u = 0, Lu = d t u + P(y, d y )u = inMxl. 

Moreover, u satisfies the estimates 

(5.4) u(x, y ,t)< ci (yoMx) V(x, t) G M x R, 

(5.5) u(x , y, t) < c 2 (x )tp*(y) V(y, t) e M x E. 

Using the semigroup property we have for all t, r > 



(5.6) / k(x,z,r)k(z,y,t + tj)dz — / k(x, z, t + tj)k(z, y, r) dz 

(5.7) = k(x, y, t + 1 + tj) 

for all j > 1. On the other hand, by (|2.13|) . the ground state 
is a positive invariant solution. Consequently, for any x G M and 
t > 0, k(x,-,r)(p G L X (M). Similarly, for all r > and y e M, 
<f*k(-,y,r) G L : (M). Hence, by estimates (I5.3P and the Lebesgue 
dominated convergence theorem, we obtain 
(5.8) 

k(x, z, t)u{z, y, t) dz = / u(x, z,t)k(z,y,r) dz = u(x,y,r + t), 
m Jm 
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where u(x,y,t) = Hindoo k(x,y,t + tj). In particular, for r = tj and 

t = —tj we have 

(5.9) 



k(x, z, tj)u(z, y, — tj) dz = / u(x, z, — tj)k(z, y, tj) dz = u(x, y, 0) 

M ' JM 

Invoking again estimates (I2.12p . and (I5.4l) -( l53j) . we see that the in- 
tegrands in (15. 9p are bounded by C(x, y)<p(z)<p*(z). Recall our as- 
sumption that ip(x)ip*(x) G L X {M), hence, the Lebesgue dominated 
convergence theorem implies 
(5.10) 



u(x, z, 0)u(z } y, 0) dz = / u(x, z, 0)u(z, y, 0) dz = u(x, y, 0), 

M JM 

where 

(5.11) u(x,y,t) := lim u(x,y,t — tj), 

j->oo 

and again we may assume that the limit in (15. lip exists. 

On the other hand, by the invariance property (I2.13p . and estimates 
fl2~T2|) . we get 
(5.12) 



u(x,z,r)(p(z)dz = cp(x), / ip*(z)u(z,y,t)dz = (p*(y). 

M JM 

In particular u(x, y, t) > 0. It follows from the Harnack inequality and 
(I5.10p that u(x,y,t) is also positive and we have 



(5.13) / / u(x, z, 0)u(z, y, 0) dzip(y) dy = / u{x,y,0)<p{y)dy. 

JM JM JM 

Consequently, (I5.12p and (15.131) imply that 



(5.14) / u(x, z,0)(p(z) dz = (p(x). 

JM 

Define integral operators 
(5.15) 



Uf(x):= / u(x,z,0)f(z)dz, U*f(y):= / u(z,y,0)f(z)dz. 

JM JM 

By (15.101) and (15.141) we see that (f(x) and u(x, y, 0) (as a function of x) 
are positive eigenfunctions of the operator U with an eigenvalue 1, and 
for every x G M, u(x,y,0) is a positive eigenfunction of U* with an 
eigenvalue 1. Moreover, for every x G M, u(x, -,0)(f G L l (M) and for 
every x,z G M, u(x, ■, 0)u(-, z, 0) G L l (M). Consequently, it follows 
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[4"8| Lemma 3.4] that 1 is a simple eigenvalue of the integral operator 
U. Hence, 

(5.16) u(x,y,0) = /3(yMx). 
By fl5TT2l) and ([5TTB1 we have 

(5.17) / (p*(z)P(y)<p(z)dz= [ ( p*(z)u(z : y 1 0)dz = cp*(y). 
Jm Jm 

Therefore, we obtain from (15. 161) and (I5.17f) that f3(y) = j v ( z )^*( z ) d z ■ 
Thus, 

(5.18) lim k{x,y,t 3 ) = u(x,y,0) = ^Y^\ A , 

3-+ 00 JmVvw (*) dz 

and the positive-critical case is proved since the limit in (I5.18P is inde- 
pendent of the sequence {tj}. 

(iii) The null-critical case: Without loss of generality, we may assume 
that PI = 0, where P is a null-critical operator in M. We need to prove 
that 

lim k^(x,y,t) = 0. 

t—^oo 

As in Lemma 14.11 (Varadhan's lemma), consider the Riemannian 
product manifold M := M x M, and let P = P Xl + P X2 be the corre- 
sponding skew product operator which is defined on M. 

If P is subcritical on M, then by part (i), lim^oo kp(x, y, t) = 0, 
and since 

kp(x,y,t) = kp(x 1 ,y 1 ,t)kp I (x 2 ,y2,t), 

it follows that lim^oo kp(x, y,t) = 0. 

Therefore, there remains to prove the theorem for the case where 
P is critical in M. Fix a nonnegative, bounded, continuous function 
/ 7^ such that <p*f G L 1 (M), and consider the solution 

v{x,t) := [ k P 1 {x,y,t)f(y)dy. 

Let t n — > oo, and consider the sequence {v n (x,t) := v(x,t + t n )}. As in 
part (ii), up to a subsequence, {v n } converges to a nonnegative solution 
ueU P {MxR). 

Invoking Lemma I4~T1 (Varadhan's lemma), we see that u(x, t) = a(t). 
Since u solves the parabolic equation Lu = 0, it follows that a(t) is a 
nonnegative constant a. 
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We claim that a = 0. Suppose to the contrary that a > 0. The 
assumption that tp*f G L l (M) and (12.1 3p imply that for any t > 

^(y)«(y, t)di/= / y*(y)(7 kp*(y,z,t)f(z)dz) dy = 
M Jm \Jm J 



(5./9)f / <p*(y)kf(y,z,t)dy)f(z)dz 
Jm \Jm J 



ip*(z)f(z) dz < oo. 



M 



On the other hand, by the null-criticality, Fatou's lemma, and (I5.19P 
we have 

oo = / (f*(z)adz= / f*{z) lim v(z, t n ) dz < 
Jm Jm n ^°° 

liminf / ip*(z)v(z, t n ) dz = I cp*(z)f(z) dz < oo. 
n ^°° Jm Jm 

Hence a = 0, and therefore 

(5.20) lim / k¥(x,y,t)f(y)dy = lim v(a;,t) = 0. 

t— >oo y w t— >oo 

Now fix y G M and let / := kp(-,y, 1). Consider the minimal solu- 
tion of the Cauchy problem with initial data /. So, by the semigroup 
property we have 

Mi 

K 

M 

h W i .. . -,-\iM i - ,. 1 i ,i - _ ;..'■/, 

M 

In view of (I2.12p (with v = 1) and (12.131) . the function / is bounded 
and satisfies f<p* G L l (M). Therefore, by (15.201) lim t _ >00 u(a;,t) = 0. 
Thus, 

lim kp(x, y, t) = lim kp(x, y, t + 1) = lim m(x, t) = 0. 

The last statement of the theorem concerning the behavior of the 
Green function G M _ x (x, y) as A — > Ao follows from the first part of the 
theorem using a classical Abelian theorem [581 Theorem 10.2]. □ 



u(x,t) := J k P (x, z,t)f(z) dz 

kp 4 (x, z, t)k% {z, y, 1) dz = k^ (x, y,t + l). 



6. Applications of Theorem 11.31 

We discuss in this section some applications of Theorem 11.31 and 
comment on related results. First we prove (11.51) . which characterizes 
Aq in terms of the large time behavior of log kp{x, y,t). 
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Corollary 6.1 ([24|). The heat kernel satisfies 

(6.1) lim \ogk¥(x,y,t) = _ Aq(P;M) Xjy€M . 

f— >-oo t 

Proof. The needed upper bound for the validity of (16.11) follows directly 
from Theorem 11.31 and (11.4j) . 

The lower bound is obtained by a standard exhaustion argument 
and Theorem 11.31 (cf. the proof of J27J Theorem 10.24]). Indeed, let 
{Mj}j S L l be an exhaustion of M. Recall that since Mj is a smooth 
bounded domain, the operator P — A (-P, Mj) is positive-critical in Mj. 
Therefore, Theorem 11.31 and the monotonicity of heat kernels with 
respect to domains (see part 4 of Lemma 12.21) imply that 

log*#(s, y ,t) > log k?(x,y,t) = _ x 

t^oo t ~ t^oo t 

Since lim^oo A (P, Mj) = Xq(P, M), we obtain the needed lower bound. 

□ 

We now use Theorem 11.31 to strengthen Lemma 14.11 More precisely, 
in the following result we obtain the large time behavior of solutions 
of the Cauchy problem with initial conditions which satisfy a certain 
(and in some sense optimal) integrability condition. 

Corollary 6.2 (|51j). Let P be an elliptic operator of the form \2. 1\) 
such that Ao > 0. Let f be a continuous function on M such that 
v*f G L 1 (M) for some v* G C P *(M). Let 

u(x,t)= [ kf{x,y,t)f(y)dy 

be the minimal solution of the Cauchy problem with initial data f on 
M. Fix K CC M. Then 



lim sup \u(x, t) — J-"(x) | = 0, 



where 



f M f(y)<p*(v) d v 



(p(x) if P is positive-critical in M, 



otherwise. 

Proof. Since v*f G L l (M), estimate (I2.12p . and the dominated conver- 
gence theorem imply that 

Urn „(,,*) = lim / *Jf(,,,, t )/ to )d,= / lim [#(,,»,*)/(»)] d», 

t^OQ t^OO J M J M t^OO 
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and the claim of the corollary follows from Theorem 11.31 □ 

Assume now that PI = and f M kp(-, y, t) dy = 1 (i.e. 1 is a 
positive invariant solution of the operator P in M). Corollary 16.21 
implies that for any j > 1 and all x G M we have 

f M I —jr 1 — — ■ — if P is positive-critical in M, 

hm / kp{x,y,t)dy=< J M ip*{y)dy 

3 I 1 otherwise. 

Suppose further that P is not positive-critical in M, and / is a bounded 
continuous function such that lim inf x ^.oo f(x) = e > 0. Then 



(6.2) lim inf / k${x, y, t)f{y) dy > e 

t— >oo 



Mi 

tv 

M 

Hence, if the integrability condition of Corollary 16.21 is not satisfied, 
then the large time behavior of the minimal solution of the Cauchy 
problem may be complicated. The following example of W. Kirsch and 
B. Simon [34] demonstrates this phenomenon. 



Example 6.3. Consider the heat equation in R d . Let R 3 := e eJ and 
let 

f(x) := 2 + if Rj < sup |^| < R j+1 , j > 1. 

l<i<d 

Let u be the minimal solution of the Cauchy problem with initial data 
/. Then for t ~ RjRj + i one has that u(0,t) ~ 2 + (— l)- 7 , and thus 
u(0, t) does not have a limit. Note that by Lemma H~Tl for d = 1, u(x, t) 
has exactly the same asymptotic behavior as u(0,t) for all x G R. 

In fact, it was proved in [56] that for the heat equation on M. d , the 
following holds: for any bounded function / defined on R , the limit 

lim / k Rd A (x,y,t)f(y)dy 

t^oo J Rd 

exists, if and only if the limit 

i im fp? 1 pm / f(y) d v 

R-^oo \B(X,R)\ J B {x,R) 

exists. Moreover, the values of the two limits are equal. For an exten- 
sion of the above result see [31J. We note that such a theorem is false 
if the average value of / is taken on solid cubes. 

The next three corollaries concern elliptic operators on manifolds. 
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Corollary 6.4 ([12]). Let M be a noncompact complete Riemannian 
manifold, and denote by A the corresponding Laplace- Beltrami opera- 
tor. Then 

(6.3) £q := lim k M A (x,y,t) = 

t— >oo 

if and only if M has an infinite volume. 

Proof. If M has a finite volume, then —A is positive-critical in M, and 
by Theorem 11.31 we have £ = (vol(M)) > 0. 

Suppose now that M has an infinite volume. If Ao > 0, Corollary 16. II 
implies that £q = 0. On the other hand, if Ao = 0, then since u — 1 is 
a positive harmonic function which is not in L 2 (M), the uniqueness of 
the ground state implies that —A is not positive-critical in M. Hence, 
by Theorem 11.31 we have £q = 0. □ 

Before studying the next two results we recall three basic notions 
from the theory of manifolds with group actions. 

Definition 6.5. Let G be a group, and suppose that G acts on M. For 
any real continuous function v and g G G, denote by v 9 the function 
defined by v 9 (x) := v(gx). Let IR* = R \ {0} be the real multiplicative 
group. 

1. A nonzero real continuous function / on M is called G -multi- 
plicative if there exists a group homomorphism 7 : G — > IR*, such that 

f(gx) = 1(9) f{x) VgeG,xeM. 

2. A G-group action on M is compactly generating if there exists 
K m M such that GK = M (see [38]). 

3. The operator P is said to be G-equivariant if 

P[u 9 } = {P[u)Y Vg e G. 

Corollary 6.6 (cf. [12J). Suppose that a noncompact manifold M is a 
covering of a compact Riemannian manifold, and consider the Laplace- 
Beltrami operator A on M . Then 

(6.4) £ := lim e Xot k^ A (x, y, t) = 0. 

t— >oo 

In particular, —A — Ao is not positive-critical in M. 

Proof. Suppose that £ > 0. By Theorem 11.31 the operator —A — Ao is 
positive-critical in M, and in particular, the ground state ip is in L 2 (M). 
The uniqueness of the ground state implies that ip is G-multiplicative, 
where G is the deck transformation. But this contradicts the assump- 
tion that p G L 2 , since G is an infinite Group. □ 
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The following corollary deals with manifolds with group actions and 
general equivariant operators, without any assumption on the growth 
of the acting groups. It generalizes Corollary 16. 6[ and seems to be new. 

Corollary 6.7. Let M be a noncompact manifold, and let G be a group. 
Suppose that G acts on M , and that the G-group action on M is com- 
pactly generating. Assume that P is a G-equivariant second-order ellip- 
tic operator of the form (12.11) which is defined on M. Suppose further 
that any nonzero G -multiplicative C 2 ' a -function is not integrable on M . 
Then 

(6.5) i := lim e Xot k^ (x,y,t) = 0. 

t— >oo 

Proof. Suppose that £ > 0. By Theorem 11.31 the operator P — A is 
positive-critical in M. In particular, dimCp-\ (M) = dim Cp*-\ (M) = 
1. By the uniqueness of the ground state, it follows that the corre- 
sponding ground states ip and ip* are G-multiplicative. Therefore, their 
product pp* is also G-multiplicative, and we arrived at a contradiction 
since by our assumption a nonzero G-multiplicative C 2,Q -function is 
not integrable on M. □ 

Remark 6.8. 1. Corollary 16.71 clearly applies to the case where M is 
a regular covering of a compact Riemannian manifold, and G is the 
corresponding deck transformation. In fact, in this case, by [38, The- 
orem 5.17], the product of the ground states is a positive G-invariant 
function, and therefore it is not integrable on M. 

2. Corollary 16 . 71 applies in particular to the case where P is defined on 
M d , and P is Z d -equivariant (that is, P has Z d -periodic coefficients). In 
this case it is known that P — Xq is (null)-critical if and only if d = 1, 2. 
For further related results, see [38] . 

Finally, let us introduce an important subclass of elliptic operators 
P such that P — Ao is positive-critical in M. 

Definition 6.9. Assume that P is symmetric and P — Ao is positive- 
critical in M with a ground state p> satisfying the normalization con- 
dition ||y||i^(Af) = 1. Let Tp (t) be the corresponding (Dirichlet) semi- 
group on L 2 (M) generated by P. We say that T^(t) is intrinsically 
ultracontr active if for any t > there exists a positive-constant Ct such 
that 

< e Xot k^(x, y, t) < ct<p{x)tp{y) Vx, y G M. 

Example 6.10. If M is a smooth bounded domain and P := —A, 
then the semigroup Tp(t) is intrinsically ultracontractive on L 2 (M). 
Also, let M := R d and P := -A + (1 + \x\ 2 ) a/2 , where a E E. Then 
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Tp(t) is intrinsically ultracontractive on L 2 (M) if and only if a > 2 
(see also [201 HI] and references therein) . 

It turns out that in the intrinsically ultracontractive case the rate 
of the convergence (as t — > oo) of the /i-transformed heat kernel is 
uniformly exponentially fast. 

Theorem 6.11 ( |20l Theorem 4.2.5]). Assume that P is symmetric 
and subcritical in M, and suppose that Tp(t) is intrinsically ultra- 
contractive on L 2 (M). Then there exists a complete orthonormal set 
{</2j}!-*L in L 2 (M) such that ipo = tp, and for any j > the function 
ifj is an eigenfunction of the Friedrichs extension of the operator P 
with eigenvalue Xj, where {Xj} is nondecreasing. Moreover, the heat 
kernel has the eigenfunction expansion (12. 9p . and there exist positive- 
constants C and 8 such that 

Mi 
»_p \ 

<Po(x)<p (y) 



c \ t kp(x, y, t) ^ 



<Ce~ st Vx,y eM,t>l. 



Furthermore, for any e > the series 

e x t kp*(x,y,t) = ^ c _ Xnt <f n (x)<f n (y) 

converges uniformly on M x M x [e, oo). 

Remark 6.12. Murata [15] proved that if Tp J (t) is intrinsically ultra- 
contractive on L 2 (M), then 1 is a small perturbation of P in M. In 
particular, for any n > the function <f n /<fo is bounded, and has a 
continuous extension up to the Martin boundary of M (with respect 
to P) [50]. On the other hand, an example of Banuelos and Davis in 
[7] gives us a finite area domain McK 2 such that 1 is a small pertur- 
bation of the Laplacian in M, but the corresponding semigroup is not 
intrinsically ultracontractive. 

Remark 6.13. In the null-recurrent case, the heat kernel may decay 
very slowly as t — > oo, and one can construct a complete Riemannian 
manifold M such that all its Riemannian products M\j > 1 are null- 
recurrent with respect to the Laplace-Beltrami operator on M- 7 (see 



7. Davies' conjecture concerning strong ratio limit 

Having proved in Section [5] that lim^oo e Xot kp (x, y, t) always exists, 
we next ask how fast this limit is approached. It is natural to conjecture 
that the limit is approached equally fast for different points x, y G M. 
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(7-1) 



Note that in the context of Markov chains, such an (individual) strong 
ratio limit property is in general not true [13]. The following conjecture 
was raised by E. B. Davies |22j in the selfadjoint case. 

Conjecture 7.1 (Davies' conjecture). Let Lu — Ut + P(x,d x )u be 
a parabolic operator which is defined on a noncompact Riemannian 
manifold M, and assume that X (P,M) > 0. Fix a reference point 
Xq G M. Then 

Xff(s,y,t) 
kp(x ,x ,t) 
exists and is positive for all x, y G M. 

The aim of the present section is to discuss Conjecture 17. H and closely 
related problems, and to obtain some results under minimal assump- 
tions. Since, the conjecture does not depend on the value of Ao, we 
assume throughout the present section that Ao = 0. 

Conjecture 17.11 was recently disproved by G. Kozma [35] in the dis- 
crete setting: 

Theorem 7.2 (Kozma [35J). There exists a connected graph G with 
bounded weights and two vertices x,y G G such that 

(7.2) hm X > n \ 

n^oo k[y,y,n) 

does not exist, where k is the heat kernel of the lazy random walk such 
that the walker, at every step, chooses with probability 1/2 to stay in 
place, and with probability 1/2 to move to one of the neighbors (with 
probability proportional to the given weights). 



In addition, Kozma indicates in [35] how the construction might be 
carried out in the category of Laplace-Beltrami operators on manifolds 
with bounded geometry. As noted in [25]) the bounded weights prop- 
erty in the graph's setting is the analogue to the bounded geometry 
property in the manifold setting, a property that in fact was not as- 
sumed in Davies' conjecture. We note that in Kozma's example the 
ratio k( x > x > n ) j s bounded between two constants that are independent of 

k(y,y,n) r 

t. 

Nevertheless, as we show below, there are many situations where 
Conjecture 17.11 holds true. 

Remark 7.3. 1. Theorem 1 1 . 31 implies that Conjecture 17. ll holds true in 
the positive-critical case, therefore, we assume throughout the present 
section that P is not positive- critical. 
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2. We also note that if P is symmetric, then the conjecture holds if 

(7.3) dimC P _ Ao (M) = 1, 

(i.e. in the "Liouvillian" case) see [U Corollary 2.7]. In particular, 
it holds true for critical selfadjoint operators. It also holds for the 
Laplace-Beltrami operator on a complete Riemannian manifold of di- 
mension d with nonnegative Ricci curvature [22J. 

3. Recently Agmon [2] obtained the exact asymptotics (in (x,y,t)) 
of the heat kernel for a Z d -periodic (non-selfadjoint) operator P on 
IR d , and for an equivariant operator P defined on abelian cocompact 
covering manifold M. In particular, it follows from Agmon's results 
that Conjecture 17.11 holds true in these cases. Note that in these cases 
(and even in the case of nilpotent cocompact covering) it is known [58] 
that 

(7.4) dimC P _ Ao (M) = dimC P *_ Ao (M) = 1. 



4. For other particular cases where the conjecture holds true see 

bq m in m [221165]. 

Remark 7.4. It would be interesting to prove Conjecture 17.11 at least 
under the assumption 

(7.5) dimC P _ Ao (M) = dimC P *_ Ao (M) = 1, 

which holds true in the critical case and in many important subcritical 
(Liouvillian) cases. For a probabilistic interpretation of Conjecture l7.lt 
see |4]. 

Remark 7.5. Let t n — > oo. By a standard parabolic argument, we 
may extract a subsequence {t nk } such that for every x,y G M and 
s < 

/ 7 m / \ i- kp(x,y,s + t n ) 

(7.6) a{x, y, s) := hm — w —— 

exists. Moreover, a(-, y, ■) e %>(MxR_), and a(x, -, •) eH P .(Mxl_). 
Note that in the selfadjoint case, we can extract a subsequence {t nfc } 
such that the limit function a satisfies a(-,y, •) G Hp(M x M) [52] . 

Remark 7.6. Consider the complete (two-dimensional) Riemannian 
manifold M that is constructed in [49]. Then M does not admit non- 
constant positive harmonic functions, Ao(— A,M) = 0. Nevertheless, 
the heat operator does not admit any Ao-invariant positive solution. In 
particular, M is stochastically incomplete (this construction disproves 
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Stroock's conjecture concerning the existence of a Ao-invariant positive 
solution). 

On the other hand, since M is Liouvillian, Remark 17.31 implies that 
Conjecture 17.11 holds true on M. Hence, the limit function 

/ n k™ A (x,y,t) 

(7 - 7) " U - !l): ^U^f,,,,,,/) 

(which equals to the constant function 1) is not a Ao-invariant posi- 
tive solution. Compare this with [221 Theorem 25] and the discussion 
therein above Lemma 26. In particular, it follows that the function 



H{x) := sup 



K p 4 (x, x, t) i<t<oQ 



(7.8) hm p ; ^;: , ;,t; ; = i v*, y e m, s g 



Kf{x Q ,x Q ,t) 

is not slowly increasing in the sense of |22j . 

Suppose that P is symmetric (and A D = 0). Using (I2.9P and a stan- 
dard exhaustion argument, it follows [22] that for a fixed x G M, the 
function t h- > kp(x,x,t) is a nonincreasing log-convex function, and 
therefore, a polarization argument implies that the following strong 
ratio property holds true. 

kp*(x,y,t + s ) 
A~ k^(x,y,t) 

In the nonsymmetric case, Corollary 16.11 and the parabolic Harnack 
inequality imply: 

Lemma 7.7 (|52j). For every x,y G M and s G K, we /iawe £/ia£ 

, 70 n ,. . f kp<(x,y,t + s) k^(x,y,t + s) 
7.9 hmmt , — < 1 < hmsup . f/ r — . 

1 ; *^oc k¥(x,y,t) ~ ~ t-+oc kp*{x,y,t) 

In particular, if lim t ^ 00 [/cp / (x, y, t + s)/kp(x, y, £)] exists, it equals to 
1. 

Remark 7.8. If there exist £o,2/o G M and < So < 1 such that 

/ 7 -,M »w x ,• kp(x , y ,t + s ) 

(7.10) M(x , ?/o, s ) := hmsup — — < oo, 

t^oo Kp [xo, yo, t) 

then by the parabolic Harnack inequality, for all x,y, z,w G K CC M, 
t > 1, we have the following Harnack inequality of elliptic type: 

(7.11) k^(z,w,t) < C 1 k^(x ,y ,t+^) < 

C 2 k P 4 (x ,yo,t-^-) <C 3 k^{x,y,t). 
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Similarly, (I7.10p implies that for all x,y G M and rGl: 

< m{x, y,r):= lim mf — < 

lim sup Wu 7, t\ = M(x ' y ' r) < 00 • 

i-»oo ftp (^0, 1/0, £J 

Note that 1 17. lOR is obviously satisfied in the symmetric case, and con- 
sequently (17. lip holds true [201 Theorem 10]. 

It turns out that the validity of the strong limit property (I7.8P (in 
the nonsymmetric case) implies the validity of Davies' conjecture if in 
addition (I7.5P is satisfied. We have: 

Lemma 7.9 Q52J). (a) The following assertions are equivalent: 

(i) For each x,y G M there exists a sequence Sj — > of negative 
numbers such that for all j > 1, and s = Sj, we have 

(7.J2) Um *iW + '> = 1. 

^ «-*» kjf(x,y,t) 



(ii) The ratio limit in \7.12) exists for any x,y G M and s6l. 
(Hi) Any limit function u(x,y,s) of the quotients k A^ ,y ' tn ~[ s } with 

fcp yXQ ,XQ )t n J 

t n — > oo does not depend on s and has the form u(x, y), where u(-, y) G 
Cp(M) for every y G M and u(x, •) G Cp*(M) for every x G M. 

(b) If one assumes further that (7T5p is satisfied, then Conjecture \7.1\ 
holds true. 

Moreover, Conjecture 7. 1 holds true if M ^ M. d is a smooth un- 
bounded domain, P and P* are (up to the boundary) smooth operators, 
(I7.12p holds true, and 

(7.13) dim C° P (M) = dim C P »(M) = 1, 

where C P (M) denotes the cone of all functions in Cp(M) that vanish 
on dM . 

Proof, (a) By Lemma [7. 7\ if the limit in (17.121) exists, then it is 1. 

(i) =r- (ii). Fix Xo, G M, and take sq < for which the limit (17.121) 
exists. It follows that any limit solution u(x,y,s) G "Hp(M x IR_) 

of a sequence with t n — > oo satisfies u(xo,y ,s + s ) = 

u(x ,y ,s) for all s < 0. So, u(x ,y ,-) is s -periodic. It follow from 
our assumption and the continuity of u that u(xq, yo, •) is the constant 
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function. Since this holds for all x,y G M and u, it follows that fIT. 12D 
holds for any x,y G M and s G K. 



ii) =>■ (iii) . Fix y G M. By Remark I7.5[ any limit function u of the 

tfoiiZ) } with *" ^ 00 belongs t 

k^(x,y,t + s) k^(x,y,t) kp*(x,y,t + s) 



sequence (tw^^tt) with t n —t 00 belongs to T-tp(M x R_). Since 

^ ftp \XQjXQ,t n ) J 



(7 ' 14) k^(x ,x Q ,t) k¥(x ,x ,t) kp*(x,y,t) ' 

(17.12p implies that such a m does not depend on s. Therefore, u 
u(x, y), where u(-, y) G Cp(M) and u(x, •) G Cp*(M). 



(7.15) 



iii) =>- (i). Write 

A;|, f (x, y,t + s) kpt(x, y,t + s) kp 1 (x , x , t) 



k^{x,y,t) k^(x ,x ,t) k A /(x,y,t) 



Let t n — > 00 be a sequence such that the sequence {^r^^j-l} 
verges to a solution in 7ip(M x R_). By our assumption, we have 

kp(x,y,t n + s) kp{x,y,t n ) . . 
hm "7m? 7T~ = lim TmT TT = w(^,J/J > °> 



COU- 



rt— >■ 00 



which together with (I7.15P implies (I7.12p for all sGl. 

(b) The uniqueness and (iii) imply that iSJr^-^; — > u } x ^Jf\ ) where 

Kp yXQ^XQjtj UyXQ )U \XQ J 

u G C P (M) and u* G C P * (M), and Conjecture O holds. □ 

Remark 7.10. Assume that one of the assumptions of part (a) of 
Lemma 17.91 is satisfied. Then by the lemma, any limit function of the 
sequence is °f the form a(x,y), where for every y G M, 

the function a(-,y) G Cp(M), and a(x, •) G Cp*(M) for every x G M. 
But in general, a(x, y) does not need to be a product of solutions of 
the equations Pu = and P*u = 0, as is demonstrated in [15], in the 
hyperbolic space, and in [521 Example 4.2]. 

In the null-critical case we have: 

Lemma 7.11 ([52]). Suppose that P is null- critical, and for each x, y G 
M there exists a sequence {sj} of negative numbers such that Sj — > 0, 
and 

(7.16) u mW MW±£)>! 



/or s = s,-, j = 1, 2, . . . . JTien Conjecture\7. 1\ holds true. 
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Proof. Let u(x,y,s) be a limit function of a sequence 

with t n — )■ oo and s < 0. By our assumption, u(x,y, s + Sj) > u(x,y, s), 
and therefore, u s (x,y,s) < for all s < 0. Thus, u(-,y,s) (respect., 
•, s)) is a positive supersolution of the equation Pit = (respect., 
P*u = 0) in M.. Since P is critical, it follows that u(-,y,s) G Cp(M) 
(respect., u(x, -,s) G Cp*(M)), and hence u s (x,y,s) = 0. By the 
uniqueness, u equals to J^yp^xo) ' anc ^ Conjecture 17.11 holds true. □ 

The large time behavior of quotients of the heat kernel is obviously 
closely related to the parabolic Martin boundary (for the parabolic 
Martin boundary theory see [TBI 121^ ^1)- The next result relates Con- 
jecture [7j] and the parabolic Martin compactification of Hp{M x 

Theorem 7.12 (|52j). Assume that (I7.10p holds true for some XQ,yo G 
M , and s > 0. Then the following assertions are equivalent: 
(i) Conjecture \7.1\ holds true for a fixed xq G M. 

(ii) 

k?(x,y,t) 



(7.17) lim 



exists, and the limit is positive for every x,y,x\,yi G M. 
(Hi) 

(7.18) lim f M [ X,V,t l and lim 

K ' t^oo kp 1 (y, y, t) hp* (x, x, t) 

exist, and these ratio limits are positive for every x,y G M. 

(iv) For any y G M there is a unique nonzero parabolic Martin bound- 
ary point y for the equation Lu = in M xl which corresponds to 
any sequence of the form {(y, —t n )} c ^ =1 such that t n — > oo, and for any 
x G M there is a unique nonzero parabolic Martin boundary point x for 
the equation u t + P*u = 0inMxR which corresponds to any sequence 
of the form {(x, — t n )}^ =1 such that t n — > oo. 



Moreover, if Conjecture 7. 1 holds true, then for any fixed y G M (re- 
spect., x G M), the limit function a(-,y) (respect., a(x, ■)) is a positive 
solution of the equation Pu = (respect., P*u = 0). Furthermore, the 



Martin functions of part (iv) are time independent, and holds 
for all x,y G M and s G JL 
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Proof, (i) =>• (ii) follows from the identity 

k*?(x,y,t) kp*{x,y,t) (k¥{x x ,y u ty 



kp(xi,yi,t) kf{x Q ,x Q ,t) V^p (x ,x Q ,t) 

(ii) =>• (iii). Take x\ = y\ = y and xi — y± — x, respectively. 

(iii) =>■ (iv). It is well known that the Martin compactification does 
not depend on the fixed reference point x$. So, fix y e M and take it 
also as a reference point. Let {— t n } be a sequence such that t n — > oo 

and such that the Martin sequence { ~fi^~7^ } converges to a Martin 

function Kp(x,y,t). By our assumption, for any t we have 

k^(x,y,t + t n ) k^(x,y,r) 
hm —TJ7 = hm -frj) r = b[x) > 0, 

n->oc (y, y, t + t n ) r^oo fc^ (y, y, r) 

where b does not depend on the sequence {— t n }. On the other hand, 



kf{y,y,t + t n ) 



Um , ; , v = ^ (yj _ t) = f{t) 

p 



»-x K%(y,y,t 



Since 



k^(x,y,t + t n ) k¥(x,y,t + t n ) fc P %y,t + t n ) 



kp*(y,y,t n ) kp*{y,y,t + t n ) k^(y,y,t n ) ' 

we have 

Kp*(x,y,t) = b(x)f(t). 

By separation of variables, there exists a constant A such that 

Pb-Xb = onM, /' + A/ = onM, /(0) = 1. 

Since 6 does not depend on the sequence {— £„}, it follows in particular, 

that A does not depend on this sequence. Thus, lim T _ ) . 00 ~far|~y — 

/(£) = e~ A *. Lemma 17.71 implies that A = 0. It follows that b is a 
positive solution of the equation Pu = 0, and 

(7.19) K»{x,y,t) = hm = b(x). 



The dual assertion can be proved similarly. 

(iv) =>- (i). Let Kp (x,y,t) be the Martin function given in (iv), and 
s > such that Kp 1 (x ,y, s /2) > 0. Consequently, Kp (x,y, s) > 
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for s > So. Using the substitution r = s + so we obtain 

fcff(z,y,r) f fcff(s,y,s + s ) 
(7.20) hm -jf- = hm <^ — — x 

kp(xo,x ,T) s^oc { kp(y,y,s) 

kp*(y, y, s) kp*(x , y, s + 2s ) \ K$f{x, y, s )K*Z (2*, y, s ) 



k¥(x ,y,s + 2so)kp(xo,x ,s + s Q ) J Kp*(x , y, 2s ) 

The last assertion of the theorem follows from (17.191) and Lemma 17.91 

□ 

Finally we mention a problem which was raised by Burdzy and Sal- 
isbury [H] for P = -A and M C R d . 

Question 7.13 Q52J). Assume that X = 0. Determine which minimal 
functions (in the sense of Martin's boundary) in Cp(M) are minimal 
in 1-Lp{M. x In particular, is it true that in the critical case, the 

ground state tp is minimal in ~Hp(Ai x M_) ? 

8. Comparing decay of critical and subcritical heat 

KERNELS 

In this section we are concerned with the large time behavior of 
the heat kernel kp with regards to the criticality versus subcriticality 
property of the operator P. Since for any fixed x, y G M, x ^ y, 
we have that kp(x,y, •) G L 1 (IR + ) if and only if P is subcritical, it is 
natural to conjecture that under some assumptions the heat kernel of 
a subcritical operator P + in M decays (in time) faster than the heat 
kernel of a critical operator Pq in M. Hence, our aim is to discuss the 
following conjecture in general settings. 

Conjecture 8.1 ([21]). Let P + and P$ be respectively subcritical and 
critical operators in M . Then 

kpi(x,y,t) 

(8.1) lim -%i ' = 

locally uniformly in M x M. 

Conjecture 18.11 was stimulated by the following conjecture of D. Kr- 
ejciffk and E. Zuazua [36J: 

Let P + and Pq be, respectively, selfadjoint subcritical and critical 
operators defined on L 2 (M,dx). Then 

(8.2) lim \\f P ^ M ' w ^ L2 ^ dx) = 

t^oo || e Fot \\L2(M,Wdx)^L 2 (M,dx) 

for some positive weight function W. 
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Remark 8.2. Theorem 11.31 implies that Conjecture EH] obviously holds 
true if Pq is positive-critical. Moreover, by Corollary 16. li the conjecture 
also holds true if X (P + ,M) > 0. Therefore, throughout this section 
we assume that Xq(P + , M) = 0, and Pq is null-critical in M. 

Example 8.3. In [UJ M. Murata obtained the exact asymptotic for 
the heat kernels of nonnegative Schrodinger operators with short-range 
(real) potentials defined on M. d , d > 1. In particular, [HI theorems 4.2 
and 4.4] imply that Conjecture 18. II holds true for such operators. 

The following theorem deals with the symmetric case. 

Theorem 8.4 Q24J). Let the subcritical operator P + and the critical 
operator Pq be symmetric in M of the form 

(8.3) Pu = -m -1 div(mAV«) + Vu. 

Assume that A + and Aq, the sections on M of End(T M) , and the 
weights m + and tuq, corresponding to P + and Pq, respectively, satisfy 
the following matrix inequality 

(8.4) m + (x)A + (x) < Cm (x)Ao(i) for a.e. x G M, 

where C is a positive constant. Assume further that for some fixed y\ G 
M there exists a positive constant C satisfying the following condition: 
for each x G M there exists T(x) > such that 

(8.5) k¥ + (x,y u t) < Ck%(x, yi ,t) Vt > T{x). 
Then 

k¥(x,y,t) 

(8.6) lim ' = 

t->oo k%(x,y,t) 

locally uniformly in M x M. 

Proof. Recall that in light of Remark 18.21 we assume that Xq(P + , M) = 
0. Suppose to the contrary that for some x ,y G M there exists a 
sequence {t n } such that t n — > oo, and 

kp + (xo,yo,t n ) 

n-X^o k^ o (x ,y ,t n ) 

Consider the sequence of functions {un}™^ defined by 



^• 7 ) J™, mrz - : x =K>0. 



kp + (x,yo,t n + s) 

i[x,s): = 

We note that 



u n (x, s) := —j^- — — x G M, s G 

kP {XQ,yo,t n ) 



_ k p+ (x, 2/o, tn + s) k^ + (a? , 2/o, tn) 
kp + {xQ,yQ,t n ) kp o (xo,yo,t n ) 
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Therefore, by assumption (18. 7p and Remark 17751 it follows that we may 
subtract a subsequence which we rename by {u n } such that 

lim u n (x, s) = u+(x, s), 

where u + G H P+ (M x R) and u + > 0. 
On the other hand, 



k^ + (x,y ,t n ) k¥ + (x,y ,t n ) 
v n [x) :— M — — ■ — u n {x, s) M ; r • 

kp Q {x ,y ,t n ) kp + (x,y ,t n + s) 

By our assumption, X (P + ,M) = and P + is symmetric, therefore 
(17.81) implies that 

kp + (x,yo,t n ) 

lim 



n->oo k^ + (x,y ,t n + s) 
Therefore, 

lim v n (x) = lim u n (x,s) = u + (x,s), 

n— >oo n— >oo 

and u + does not depend on s, and hence u + is a positive solution of 
the elliptic equation P + u = in M and we have 

kp (x,y ,t n ) 
(8.8) lim —TjT r = u+(x). 

On the other hand, by Remark 17.31 we have 



kp {x,yp,t n ) <p(x) 



where <p is the ground state of Pq. 
Combining (I8.8p and (18.91) . we obtain 

k p (^yo^) | uo(x) 

On the other hand, by assumption (18. 5p and the parabolic Harnack 
inequality there exists a positive constant C% which depends on P + , 
P , yo, and y\ such that 

(8.11) C^ + {x,y Q ,t-l) <k^ + {x,y u t) 

< Ck%(x, yi ,t) < CC 1 k%(x,y Q ,t + 1) VxeM,t> T(x). 
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Moreover, by f!7.8p we have 
(8.12) 



kf + {x,y ,t-l) k%(x,y ,t + l) 
tj—, ; — = 1, and hm — hn-, ; — 



lim — pjrjr- — = 1 and lim — ^tt s — = 1 Vx G M. 

^co k¥ + (x,y ,t) t^oo k%(x,y ,t) 



Therefore, (18.111) and (18. 12)) imply that there exists C > such that 

(8.13) k¥ + {x,y ,t) < C k%(x,y ,t) Vx G M, t > T{x). 

Consequently, (18. ldj) and (18.131) imply that 

u +{ x ) < CqUo(x) = C ip(x) Wx G M. 

Therefore, using (18.41) we obtain 
(8.14) 

(w+) 2 (x)m + (x)A + (x) < C2(p 2 {x)m (x)A (x) for a.e. x G M, 

where C2 > is a positive constant. Thus, Theorem 12.111 implies that 
P + is critical in M which is a contradiction. The last statement of the 
theorem follows from the parabolic Harnack inequality and parabolic 
regularity. □ 



By the generalized maximum principle, assumption (18. 5ft in Theo- 
rem 18.41 is satisfied with C = 1 if P + = Pq + V, where Pq is a critical 
operator on M and V is any nonnegative potential. Note that if the 
potential is in addition nontrivial, then P + is indeed subcritical in M. 
Therefore, we have 

Corollary 8.5 (|24j). Let Pq be a symmetric operator of the form (18.31) 
which is critical in M , and let P + := Pq + V\, where V\ is a nonzero 
nonnegative potential. Then 

k¥(x,y,t) 

8.15 lim ,^t; ,y ' [ =0 

t->™ k™(x,y,t) 

locally uniformly in M x M. 

Next, we discuss the nonsymmetric case. We study two cases where 
Davies' conjecture implies Conjecture 18.11 First, we show that in the 
nonsymmetric case, the result of Corollary 18.51 for positive perturba- 
tions of a critical operator Pq still holds provided that the validity of 
Davies' conjecture (Conjecture 17. II) is assumed instead of the symmetry 
hypothesis. More precisely, we have 

Theorem 8.6 Q24J). Let Pq be a critical operator in M , and let P + = 

P + V, where V is any nonzero nonnegative potential on M. Assume 



38 YEHUDA PINCHOVER 

M 



that Davies ' conjecture ( Conjecture \ 7.1\ ) holds true for both kp and 

k?Ax,y,t) n 



k^ + . Then 



(8.16) lim , . , - w 

locally uniformly in M x M. 

Proof. Recall that in light of Remark 18.21 we assume that Xq(P + , M) = 
0. Suppose to the contrary that for some x ,yo G M there exists a 
sequence {t n } such that t n — > oo and 

kp, (x , y , t n ) 
8.17 lim ' U0 ' n} =K>0. 

Consider the functions v+ and vq defined by 

kP + {x ,y ,t) kp (x ,y ,t) 
By our assumption, 

lim v + (x,t) = u + (x), lim v (x,t) = u (x), 

t—>oo t— >oo 

where u + e C P+ (M) and u £ C Po (M). 

On the other hand, by the generalized maximum principle 

k¥(x,y ,t) 



Therefore, 



k^ + (x,y ,t n ) 

kp + (x ,y ,t n ) k¥ + (x ,yo,t n ) _ kp* + (x,y ,t n ) 
Letting n — > oo we obtain 

(8.20) Ku + (x) < u (x) xeM. 

It follows that v(x) := Uq(x) — Ku + (x) is a nonnegative supersolution 
of the equation Pqu = in M which is not a solution. In particular, 
v ^ 0. By the strong maximum principle v(x) is a strictly positive 
supersolution of the equation P u = in M which is not a solution. 
This contradicts the criticality of Pq in M. □ 

The second nonsymmetric result concerns semismall perturbations. 
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Theorem 8.7 (|24|). Let P + and Pq = P + + V be a subcritical oper- 
ator and a critical operator in M, respectively. Suppose that V is a 
semismall perturbation of the operator PI in M. Assume further that 
Davies' conjecture ( Conjecture 7.1 ) holds true for both kp and kp 1 and 
that (18. 5ft holds true. Then 



(8.21) lim = 

e Po ( 

locally uniformly in M x M. 



kX(x,y,t) 



Proof. Recall that we (may) assume that Xq(P + ,M) = 0. Assume to 
the contrary that for some xq, yo G M there exists a sequence {t n } such 
that t n oo and 

kp (Xn, Vn, t n ) 

(8.22) lim n > =K>0. 

n-^oo kpl(x ,y ,t n ) 

Consider the functions v + and v defined by 
(8.23) 

k^ + (x,y ,t) k%(x,y ,t) c w_ n 

kp + {x ,yo,t) kpl{x ,yo,t) 
By our assumption, 

lim v + (x,t) = u + (x), lim v (x,t) = u (x), 

t—too t— »oo 

where u + G C P+ (M) and u G C Po (M). 

On the other hand, by assumption (18.51) we have for t > T(x) 
(8.24) 



kg + (a>,yo,t) 



k^ + (x,y ,t) _ k¥ + (x,y u t) g^ggj fc%Qr,a/o,t) fcg(x,yx,t) 
kpl(x,y ,t) ~ kpl(x, yi ,t) X " iM) 

kP Q (x,y u t) 

By our assumption on Davies' conjecture, we have for a fixed x 
rao _, .. *ff>,2/o,*) ^(s/p) kpl(x, yi ,t) u*( yi ) 



t^ook^ + {x,y X) t) u* + (yi) ' t-^oo kpl(x,y ,t) u^(y o y 

where and «g are positive solutions of the equation P£« = and 
P *m = in M, respectively. By the elliptic Harnack inequality there 
exists a positive constant C\ (depending on Pf, P * ,yo,y± but not on 
x) such that 

(8.26) < Ci, #4 < C X . 
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Therefore, (ESQ) and <K26h imply that 

kp (x, ynA n ) 

8.27 ,^t; [ < 2CC\ 

for n sufficiently large (which might depend on x). 
Therefore, 

( v k^ + (x ,y ,t n ) kp + (x Q , yo ,t n ) k¥ + {x,y ,t n ) 

[ } k%(xo,yo,t n ) X ^jgAJ " k%(x,y ,t n ) ~ • 

Letting n — > oo and using ( I8.22p and (18. 23 p . we obtain 

(8.29) < 2CC 1 2 m (x) xGM. 

On the other hand, since V is a semismall perturbation of P£ in M, 
Theorem 12.131 implies that uq(x) x Gp + (x,yo) in M \ B(y ,5), with 
some positive 5. Consequently, 



(8.30) u+(x) < C 2 Gf + (x,y ) x eM\B(y ,5) 

for some C 2 > 0. In other words, u + is a global positive solution of the 
equation P + u = in M which has minimal growth in a neighborhood of 
infinity in M. Therefore u + is a ground state of the equation P + u = 
in M, but this contradicts the subcriticality of P + in M. □ 



9. On the equivalence of heat kernels 

In this section we study a general question concerning the equiva- 
lence of heat kernels which in turn will give sufficient conditions for 
the validity of the boundedness condition ( 18 .5p that is assumed in the- 
orems 18.41 and 18.71 

Definition 9.1. Let Pi, i — 1,2, be two elliptic operators of the form 
(EH) that are defined on M and satisfy A (-Pj, M) = for i = 1, 2. 
We say that the heat kernels kp x (x,y,t) and kp 2 (x,y,t) are equivalent 
(respectively, semiequivalent) if 

kp x -k% on M x M x (0, oo) 

(respectively, 

fc£(-, 2/o, ■)**#(•, j/0,0 on Mx(0,oo) 



for some fixed y e M). 
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There is an intensive literature dealing with (almost optimal) condi- 
tions under which two positive (minimal) Green functions are equiva- 
lent or semiequivalent (see [3j H3j HHJ [50] and the references therein). 
On the other hand, sufficient conditions for the equivalence of heat 
kernels are known only in a few cases (see [39J, 0OJ [67]). In particular, 
it seems that the answer to the following conjecture is not known. 

Conjecture 9.2 ([24]). Let Pi and Pi be two subcritical operators of 
the form ( 12.1 p that are defined on a Riemannian manifold M such that 
Pi = P2 outside a compact set in M , and \o(Pi, M) = for i = 1, 2 . 
Then kp^ and kp 2 are equivalent. 

Remark 9.3. Suppose that Pi and P2 satisfy the assumption of Con- 
jecture I9T21 B. Devyver and the author proved recently that there exists 
C > such that 



(9.1) 

C- x G%_ x {x,y) < G%_ x (x,y) < CG%_ x (x,y) Wx, y G M, and A < 0. 



Clearly, by (11. TJ , the above estimate (19.11) is a necessary condition for 
the validity of Conjecture 19.21 

It is well known that certain 3-G inequalities imply the equivalence 
of Green functions, and the notions of small and semismall perturba- 
tions is based on this fact. Moreover, small (respectively, semismall) 
perturbations are sufficient conditions and in some sense also neces- 
sary conditions for the equivalence (respectively, semiequivalence) of 
the Green functions [HI H2J |5D]. Therefore, it is natural to introduce 
an analog definition for heat kernels (cf. [67]). 

Definition 9.4. Let P be a subcritical operator in M. We say that 
a potential V is a k-bounded perturbation (respectively, k-semibounded 
perturbation) with respect to the heat kernel kp(x,y,t) if there exists 
a positive constant C such that the following 3-k inequality is satisfied: 



for all x,y G M (respectively, for a fixed y G M, and all x G M) and 



The following result shows that the notion of /c-(semi)boundedness 
is naturally related to the (semi)equivalence of heat kernels. 

Theorem 9.5 ([24]). Let P be a subcritical operator in M, and as- 
sume that the potential V is k-bounded perturbation (respectively, k- 
semibounded perturbation) with respect to the heat kernel kp(x,y,t). 
Then there exists c > such that for any \e\ < c the heat kernels 



(9.2) 




t > 0. 
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kp +eV (x,y,t) and kp (x,y,t) are equivalent (respectively, semiequiva- 
lent). 

Proof. Consider the iterated kernels 

k?(x,y,t) j = 0, 

Jo Im k p' 1] (^ *, t - s)V(z)k¥(z, y, s) dzds j > 1. 
Using ( 19. 2ft and an induction argument, it follows that 



k$\x,y,t) :-- 



J2W\^\x,y,t)\ 

3=0 

< (1 + C\e\ + C 2 \e\ 2 + . . . ) y, t) = y— ^ k^(x, y, t) 

provided that |e| < C^ 1 . Consequently, for such e the Neumann series 

oo 
3=0 

converges locally uniformly in M x M x K + to fcp f | _ el/ (x, ?/, t) which in 
turn implies that Duhamel's formula 

(9.3) k^ +eV (x,y,t)=k^(x,y,t)- 

[ [ k P d (x,z,t-s)V(z)k P d +£V (z,y,s)dzds 

J0 J M 



e 

is valid. Moreover, we have 

kp* +£V (x,y,t) < 1 _ 1 C >, k¥(x,y,t) 
The lower bound 

C x kp*{x,y,t)<kp* +£V (x,y,t) 

(for |e| small enough) follows from the upper bound using (19.31) and 
(P- □ 

Corollary 9.6 ([24J). Assume that P and V satisfy the conditions of 
Theorem \9.5[ and suppose further that V is nonnegative. Then there 
exists c > such that for any e > — c the heat kernels kp +eV (x,y,t) 
and kp \x,y,t) are equivalent (respectively, semiequivalent). 

Proof. By Theorem 19.51 the heat kernels kp +sV (x,y,t) and kp(x,y,t) 
are equivalent (respectively, semiequivalent) for any \e\ < c. 
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Recall that by the generalized maximum principle, 

(9.4) kp* +EV (x,y,t)<kp*(x,y,t) Ve > 0. 

On the other hand, given a potential W (not necessarily of definite 
sign), and < a < 1, denote P a := P + aW, and assume that Pj > 
in M for j = 0, 1. Then for < a < 1 we have P a > in M, and the 
following inequality holds |47] 
(9.5) 

k¥ a (x,y,t) < [k^y^-^ix^t)] Vx,y e M, t > 0. 

Using (19.41) and (19. 5p we obtain the desired equivalence kp +eV x kp 
also for e > c. □ 

Finally we have: 

Theorem 9.7 Q24J). Let Po be a critical operator in M. Assume that 
V = V + — y_ is a potential such that V± > and P + := P + V is 
subcritical in M. 

Assume further that V_ is k-semibounded perturbation with respect to 
the heat kernel kp + (x, yi, t). Then condition (18.51) is satisfied uniformly 
in x. That is, there exist positive constants C and T such that 

(9.6) k^ + {x, yi ,t)<Ck%{x, yi ,t) VxGM,t>T. 

Proof. It follows from Corollary 19.61 that the heat kernels kp (x,yi,t) 
and. kp^_^_y (x, y x ,t) are semiequivalent. Note that P + + V- = P + V + , 
Therefore we have 

C- 1 k P 4 + (x,y 1 ,t)<kp* ++v _(x,y 1 ,t) 

= k% +v+ (x, yi ,t) < kpl(x, Vl ,t) VxeM,t>0. 

□ 
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